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7. A 


OGRESS  IN  TECHNICAL  APPROACH 


New  techniques  for  spectroscopy  in  supersonic  expansions 
and  in  the  development  of  specific  supersonic  sources  were 
developed. , 


\/ 

Vacuum  Ultraviolet  Absorption  Spectroscopy  in  Supersonic 
Expansions.  We  have  combined  the  techniques  of  vacuum  ultra- 
vToTetT  (VUV)  spectroscopy  together  with  planar  supersonic 
jets,  which  allows  for  the  interrogation  of  absorption  spectra 
of  large  molecules  cooled  in  supersonic  expansions  in  the  near 
VUV  region. \  The  experimental  setup  consists  of  a  high  pres¬ 
sure  Xe  lamp,  CaF2  optics,  a  vacuum  ultraviolet  spectrograph 
and  a  nozzle  slit  (0.27x90  mm,  repetition  rate  9  Hz  and  gas 
pulse  dura/ion  300  psec).  The  characteristics  of  this  spectro¬ 
scopic  sptup  are:  (i)  Energy  range  6-10  eV.  (ii)  Spectral 
resolution  0.1  A.  (iii)  Routine  measurements  of  high-energy 
absorption  spectra,  (iv)  Interrogation  of  fluorescence  excita¬ 
tion^  spectra  of  the  parent  molecule  or  its  photoproducts  with 
limiting  quantum  yields  of  Y  z  10"4. 

4/ 

Development  of  Conical  Nozzles  ;£or  Supersonic  Jets.  Conical 
nozzles  (nozzle  opening  angle  (§)=■  30^,  and  nozzle  di ameter 
D  =  0.3  mm)  were  constructed  and  used  in  conjunction  with  a 
magnetic  pulsed  valve.  The  use  of  conical  nozzles  considerably 
enhances  clustering  in  supersonic  expansions,  facilitating 
studies  of  large  van  der  Waals  complexes  and  clusters. 


ACCOMPLISHMENTS  OF  OBJECTIVES 


Energetics  of  Rydberg  States  of  Jet  Cooled  Molecules.  VUV 
absorptionspectra  of  benzene,  benzene-06  and  naphthalene 
cooled  in  planar  supersonic  expansions  were  measured  over  the 
range  2000-1600  A,  providing  evidence  on  energetics,  line — 


broadening  and  interference  effects. 
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Rydberg  States  of  Anthracene.  Several  Rydberg  transitions 
of  jet-cooled  anthracene  were  observed  in  the  spectral  region 
1900-2000  A.  These  Rydbergs  are  superimposed  on  a  broad 
"background"  of  mr*  transitions.  Up  to  now  no  fluorescence 
from  Rydberg  states  of  polyatomic  molecules  has  been  reported. 

The  1997  A  n=3  Rydberg  of  anthracene  reveals  fluorescence 
with  a  quantum  yield  of  ^  5%.  This  fluorescence  originates 
from  internal  conversion  R  vw  Sx  followed  by  Sx  — ►  S0  emission, 
as  documented  by  the  spectral  distribution  of  the  dispersed 
fluorescence. 


Intramolecular  Relaxation  of  Rydberg  States.  Information  on 
intramolecular  dynamics  of  extra valence  excitations  of  benzene 
was  obtained  from  lineshape  analysis.  The  lineshape  of  the 
3Pxy(°)  Rydberg  is  Lorentzian,  whose  homogeneous  width  result 
in  the  lifetime  t  =  0.19±0.02  psec  for  C^Hg  and  t  =  0.22±0.02 
psec  for  C6D6.  Semiquantitative  information  on  the  lifetimes 
of  some  Rydbergs  of  anthracene  was  obtained,  which  fall  in  the 
range  of  ^  0.1  psec.  These  lifetime  data  imply  that  (i)  the 
relaxation  of  the  Rydberg  is  characterized  by  moderate  energy 
gap(s),  and  Cu)  the  electronic  relaxation  rate  of  the  Rydberg 
is  considerably  less  efficient  than  that  of  the  intravalence 
excitation  in  the  same  energy  domain. 


Interference  Effects  between  Extravalence  and  Intravalence 
Molecular  Excitations.  We  have  searched  for  Rydberg-valence 
interference  effects  in  large  molecules,  which  are  expected  to 
be  exhibited  in  asymmetric  Fano-type  lineshapes  in  absorption. 
For  "isolated"  Rydbergs,  which  are  superimposed  on  a  inr*  tran¬ 
sition  in  benzene,  naphthalene  and  anthracene,  the  absorption 
lineshapes  are  symmetric  Lorentzians.  No  line  asymmetry  and 
no  antiresonances  characteristic  of  Fano  profiles  were  found 
in  that  case.  The  absence  of  interference  effects  in  this  case 
reveals  the  manifestation  of  random  interstate  coupling  of  R 
with  the  it*  manifold.  The  random  coupling  erodes  all  inter¬ 
ference  effects.  The  situation  is  different  for  nearly  lying 
Rydbergs,  where  pronounced  R-R-valence  interference  effects 
were  observed  in  the  absorption  spectra  of  jet-cooled  naphtha¬ 
lene  in  the  spectral  region  1600-1650  A,  providing  information 
on  the  homogeneous  contribution  to  high-energy  molecular  coupl¬ 
ing  phenomena . 


Rotational  State  Dependence  of  Intramolecular  Dynamics.  Rota¬ 
tional  effects  on  interstate  coupling  are  of  considerable  cur¬ 
rent  interest.  Absolute  fluorescence  quantum  yields  from  photo- 
selected  rotational  states  were  measured  for  the  electronic 
origins  of  the  Sj  state  of  pyrazine.  Strong  rotational  state 
dependence  was  observed  providing  novel  information  on  inter¬ 
state  coupling  for  the  intermediate  level  structure.  This 
unique  information  cannot  be  extracted  from  time-resolved  decay 
lifetimes. 

The  Coupling  between  Intrastate  Vibrational  Energy  Redistribution 
and  Interstate  Electronic  Relaxation.  We  have  documented  some 
universal  characteristics  of  the  decay  lifetimes  and  fluorescence 
quantum  yields  from  the  S2  manifold  of  large  molecules,  which 
originate  from  the  coupling  between  intrastate  vibrational  energy 
redistribution  and  interstate  electronic  relaxation.  The  time- 
resolved  total  fluorescence  decay  excited  by  a  psec  laser  from  the 
Si  state  of  jet-cooled  9CN-anthracene  exhibits  nonexponential  decay 
in  the  energy  range  Ey  =  1200-1740  cm"1  above  the  Sx  origin,  which 
does  not  originate  from  dephasing  but  rather  manifests  the  effects 
of  intrastate  intermediate  level  structure  for  vibrational  energy 
redistribution  on  intersystem  crossing. 


Fluorescence  Quantum  Yields  for  Highly-Excited  States  of  Large 
Molecules'!  Fluorescence  quantum  yields  Y  from  high  electronic 
excitations  of  naphthalene,  anthracene  and  tetracene  in  the  energy 
range  5. 0-6. 5  eV  have  been  recorded.  Y  exhibits  an  exponential 
dependence  on  the  excess  vibrational  energy  with  the  slope  decreas¬ 
ing  with  increasing  size  of  the  molecule,  i.e.,  the  vibrational 
density  of  states.  These  data  are  characteristic  of  internal  con¬ 
version  from  Sj  to  S0.  These  results  are  of  interest  regarding 
recent  astrophysical  implications  of  high-energy  photophysics  of 
large  aromatic  hydrocarbons  in  outer  space. 


Photoisomerization  Dynamics  of  Trans-Stilbene  and  of  Cis-Stilbene. 
Time-resolved  fluorescence  lifetimes  from  photoselected  states  of 
trans-stilbene  were  recorded  by  the  techniques  of  picosecond  spec¬ 
troscopy  in  jets  using  a  mode-locked  dye  laser  and  a  fast  photon 
counting  system.  Decay  lifetimes  as  short  as  100±30  psec  were 
recorded.  Extensive  information  on  the  energy  dependence  of  the 


isomerization  rates  of  alkyl  stilbenes  was  obtained,  providing 
information  on  the  role  of  intramolecular  vibrational  distribu¬ 
tion  on  the  photochemistry  in  an  isolated  molecule.  The  absorp¬ 
tion  spectrum  of  jet-cooled  cis-stilbene  is  broad  and  no  vibra¬ 
tional  structure  could  be  resolved  due  to  the  congestion  of 
broadened  low-frequency  vibrational  excitation.  From  the  low 
quantum  yield  Y  $  2xl(T4  we  infer  a  Sx  lifetime  of  t  j  0.4 
psec,  which  indicates  ultrafast  relaxation  of  this  molecule. 


7. I  Photoisomerization  Dynamics  of  Alkyl  Substituted  Trans-Sti lbene. 
A  central  issue  pertaining  to  isolated-molecule  photochemistry 
involves  the  role  of  intramolecular  vibrational  relaxation  in 
determining  intramolecular  dynamics.  We  have  attempted  to  in¬ 
crease  the  density  of  vibrational  states  by  alkyl  substitution 
of  trans-stilbene  and  have  explored  the  isomerization  dynamics 
by  picosecond  time-resolved  spectroscopy.  We  have  obtained  the 
counter-intuitive  result  that  alkyl  substitution  of  trans- 
sti lbene  enhances  the  photoisomerization  rates,  while  general 
arguments  based  on  the  role  of  IVR  and  the  implications  of 
statistical  theories  indicate  that  the  rate  should  be  retarded. 
The  results  can  be  accounted  for  by  the  modification  of  the 
molecular  parameters,  i.e.,  threshold  energy,  by  alkyl  substi¬ 
tution. 


7.J  Energy-Resolved  Photoisomerization  Rates.  The  dynamics  of  the  Sj 
and  S2  electronically  excited  singlet  states  of  diphenylbutadiene 
was  interrogated  by  fluorescence  quantum  yield  measurements  over 
the  very  broad  energy  domain  of  0-7500  cm-1  above  the  (false) 
origin.  The  issues  of  the  lack  of  mode  selectivity  and  the  appli 
cability  of  statistical  theories  for  the  description  of  isolated- 
molecule  photochemistry  were  explored. 


7.K  van  der  Waals  Complexes  of  Porphyrins.  Excited-state  energetics 
and  dynamics  of  large  complexes  consisting  of  porphyrins  bound  to 
rare-gas  atoms  were  explored.  Detailed  spectroscopic  information 
on  the  S0  -*■  and  S0  -*■  S2  transitions  of  the  free-base  porphine- 
Ar  complex  was  obtained,  providing  insight  into  the  structure  of 
this  complex  as  well  as  resulting  in  a  novel  mechanism  for  micro¬ 
scopic  solvent  shifts  induced  by  configurational  distortions, 
which  are  due  to  complexing. 
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7.L  Coupling  between  Intramolecular  and  Intermolecular  Nuclear  Motion 
in  Complexes.  Intermolecular  vibrations  of  large  van  der  Waals 
complexes,  which  involve  the  motion  of  the  ligand  relative  to  the 
large  molecule,  provide  an  analogue  for  surface  vibrational  motion 
in  a  finite  system  and  constitute  the  precursors  of  phonon  modes 
in  condensed  phases.  Information  on  the  coupling  between  inter¬ 
molecular  and  intramolecular  vibrational  motion  was  obtained  for 
the  trans-stilbene*Ar  complex. 


7.M  Electron  Localization  in  Clusters.  Small  clusters  exhibit  unique 
physical  and  chemical  phenomena,  which  are  both  of  fundamental  and 
technological  significance,  and  provide  ways  and  means  to  explore 
the  "transition"  from  molecular  to  condensed-matter  systems.  We 
have  provided  a  theoretical  study  of  the  structure,  energetics  and 
dynamics  of  an  excess  electron  interacting  with  an  alkali-halide 
cluster,  which  was  explored  by  the  quantum  path  integral  molecular 
dynamics  method.  These  studies  establish  various  compositional, 
structural  and  size  dependence  of  bulk  and  surface  localization 
mechanisms  of  the  dynamic  process  induced  by  electron  attachment. 


7.N  Vibrational  Predissociation  Induced  by  Exciton  Trapping  in  Rare-Gas 
Clusters.  The  dynamics  of  exciton  trapping,  vibrational  energy 
transfer  and  vibrational  predissociation  in  an  electronically  excited 
state  of  Ar13  clusters  was  explored  by  classical  molecular  dynamics. 
This  study  constitutes  an  application  of  this  technique  for  the 
dynamics  of  electronically  excited  states  of  large  systems.  New 
mechanisms  of  ultrafast  10  psec)  vibrational  energy  flow  induced 
by  short-range  repulsion  were  documented.  In  these  systems  vibra¬ 
tional  energy  redistribution  does  not  occur  and  mode  selective  excita¬ 
tion  prevails.  In  small  (n=13)  clusters  the  consequences  of  vibra¬ 
tional  energy  flow  results  in  reactive  vibrational  predissociation, 
leading  to  the  "evaporation"  of  A r  atoms  for  the  cluster.  In  layer 
(n-55)  clusters  a  transition  from  molecular- type  reactive  behavior 
to  nonreactive  vibrational  relaxation,  which  is  characteristic  of 
condensed  phases,  was  exhibited. 
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Vacuum  ultraviolet  absorption  spectroscopy  in  supersonic  expansions 
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Remarkable  progess  has  been  accomplished  in  the  area 
of  molecular  spectroscopy  in  supersonic  expansions,  which 
unveiled  many  facets  of  excited-state  energetics  and  dynam¬ 
ics  of  large  molecules  excited  in  the  energy  range  0.2-6 
eV.  1~‘  Photoselective  molecular  excitation  in  jets  above  6  eV 
is  expected  to  open  new  horizons  in  the  areas  of  extravalence 
spectroscopy,  nonreactive  and  reactive  intramolecular  dy¬ 
namics,  and  state  selective  molecular  photochemistry.5  A 
significant  step  in  this  direction  was  undertaken  recently  by 
three-  and  four-photon  ionization  studies6-9  in  jets  and  by 
femtosecond  time-resolved  investigations10  of  Rydberg 
states.  We  wish  to  report  on  the  merger  between  the  tech¬ 
niques  of  vacuum  ultraviolet  (VUV)  spectroscopy  and  super¬ 
sonic  jets,  which  allows  for  the  interrogation  of  absorption 
spectra  of  large  molecules  in  supersonic  expansions  in  the 
near  VUV  (6-9  eV)  region.  The  techniques  of  absorption 
spectroscopy  in  seeded,  planar,  pulsed,  supersonic  expan¬ 
sions"-13  were  adopted  for  VUV  absorption  studies.  Planar 
supersonic  jets  of  Ar  seeded  with  benzene,  benzene-D6  or 
naphthalene  were  expanded  from  a  nozzle  slit  (dimensions 
0.27x90  mm,  repetition  rate  9  Hz  and  gas  pulse  duration 
300  p.%)  and  the  characteristic  rotational  temperatures  were 
Tr  ~20  K.  Light  (spectral  range  1 500-2000  A)  from  a  modi¬ 
fied  short-arc  high-pressure  Xe  lamp  (EG  &  G  Model  FX 
193U1,  which  was  fitted  with  a  sapphire  or  CaF;  window, 
was  focused  with  an/=  4  CaF;  lens  onto  the  jet  parallel  to 
the  nozzle  slit  at  a  distance  x  =  8  mm.  The  light  beam  was 
split  by  a  CaF:  plate  and  monitored  by  two  photomultipliers 
coated  with  a  (tetraphenyl  butadiene)  convector,  both  before 
and  after  crossing  the  planar  jet.  The  absorption  AI  /I0  was 
normalized  to  the  incident  light  intensity. 

VUV  absorption  spectra  of  jet-cooled  aromatic  mole¬ 
cules  provided  information  on  intramolecular  dynamics  of 
odd-parity  Rydberg  states  and  an  interference  between  Ryd¬ 
berg  and  intravalence  excitations.  The  electronic  origin  of 
the  3  Pxy  ['A  iu )  Rydberg  excitation  [$PX  y  (0)]  of  benzene  and 


of  benzene-£>6  (Fig.  1)  is  superimposed  on  the  XA  ,g  —  *£  lu 
intravalence  transition.1415  The  energy  of  3£[>(0)  (55  877 
cm-1  for  C(,H6and  55  946  cm-1  for  C6D6)  is  in  accord  with 
one-1415  and  three-photon16  bulb  spectroscopy.  The  line 
shape  of  3£*y(0)  is  a  Lorentzian  (Fig.  1)  with  a  width 
(FWHM)  of  A  =  32.7  +  3.1  cm-1  for  C6H6  and 
A  =29.0  +  2.0  cm-1  for  C6D6.  The  homogeneous  line- 
widths  r  and  lifetimes  r  obtained  after  correction  for  instru¬ 
mental  resolution  are  r  =  28  +  3  cm  - 1  (r  =  0. 1 9  ±  0.02  ps) 
for  C6H6,  and  £  =  24  +  2  cm-1  (r  =  0.22  +  0.02ps)  for 
C6D6.  Our  value pf  r  for  3 Pxy  (0)  is  close  to  the  linewidths  of 
27-33  cm  - 1  reported6  for  the  vibronically  induced  two-pho¬ 
ton  transitions  of  this  Rydberg  transition.  A  small  (  —  10%) 
normal  deuterium  isotope  effect  on  the  relaxation  rate  may 
exist,  which  is,  however,  within  our  experimental  uncertain¬ 
ty,  implying  that  the  electronic  relaxation  process  of  3 Pxy  (0) 
is  characterized  by  moderate  energy  gap(s).  The  homogen¬ 
eous  line  broadening  of  the  'A  lg  — <•  ‘£lu  intra  valence  excita¬ 
tion  of  jet-cooled  benzene  exceeds  £  for  3  Pxy  (0)  by  about  one 
to  two  orders  of  magnitude,  which  is  in  accord  with  old  argu¬ 
ments,  1 7  so  that  the  intramolecular  relaxation  of  the  intra  va¬ 
lence  1 E , „  excitation  presumably  occurs  within  —  10  fs,  i.e., 
on  a  time  scale  comparable  with  a  vibrational  period.  The 
Lorentzian  line  shape  of  'A  lg  — *  iPx  v  (0)  provides  compell¬ 
ing  evidence  for  the  absence  of  Rydberg-intravalence  inter¬ 
ference  effects 1 5' 1  1 9  in  benzene,  which  is  in  accord  with  Rice 
et  a/.15  The  absence  of  Fano  interference  effects20  can  be 
rationalized  within  the  framework  of  a  random  coupling 
model21 22  for  the  3£x>,(0)-1£|1,  interstate  coupling  terms 
and  for  the  'A  lg  — <•  '£,„  transition  moments  to  the  homo¬ 
geneously  broadened  background  manifold,  which  erode  all 
interference  effects.2122  Pronounced  Rydberg-intravalence 
interference  effects18  are  exhibited  in  the  spectrum  of  jet- 
cooled  naphthalene  (Fig.  2).  The  excessive  line  broadening  of 
the  resonance,  e.g.,  £=80  cm-1(r=70  fs)  for  the 
R  5  +  1360  transition  (Fig.  2),  is  close  to  the  corresponding 


FIG.  1.  Absorption  spectra  of  the 
'A  ,,  — -  iPxy[A  j„  |  electronic  origin  of 
benzene  and  benzene-£>ft  in  pulsed  planar 
supersonic  expansions  of  Ar  at  p  =  105 
Torr.  The  spectral  resolution  is  0.45  A. 
The  base  line  is  marked  as  zero  on  the  in¬ 
tensity  scale.  The  solid  curves  represent 
fits  of  the  line  shape  L  W  )  to  a  Lorentzian 
superimposed  on  a  constant  background 
L  K  |  =  A  +  C/[\A.  -  a0|-  +  (J  /2|-'].  For 
C„H„:  A  —  5.7,  C=I7.2  (in  arbitrary 
unitsl,  /i„  =  1789.65  A  and  A  —  1.08  A 
(33.7  cm  " '|  while  for  C*D,,:  A  =  6.5, 
C=  1 2.6  (in  arbitrary  unitsl.  A.,,  =  1787  45 
A  and  J  =  0.88  A  127.5  cm  '). 
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FIG.  2.  The  absorption  spectra  of  naphthalene  in  pulsed  planar  expansions 
of  Ar  ( p  =  105  Torr)  in  the  region  1630-1650  A.  The  spectral  resolution  is 
0.5  A.  The  antiresonances  corresponding  to  the  R  ’  +  1360  and  to  the  R  ‘ 
Rydberg  states  (according  to  the  assignment  of  Refs.  15and  19)  are  marked 
by  arrows. 

room  temperature  value  r=15  cm'1,15'1'1  whereupon  the 
excessive  homogeneous  broadening  make  the  thermal  con¬ 
gestion  effects.  Our  spectroscopic  approach  can  readily  be 
extended  throughout  the  LiF  region  (up  to  11  eV)  and  will 
contribute  towards  the  elucidation  of  the  microscopic 
aspects  of  intramolecular  dynamics  and  photochemistry. 
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Remarkable  progess  has  been  accomplished  in  the  area 
of  molecular  spectroscopy  in  supersonic  expansions,  which 
unv  eiled  many  facets  of  excited-state  energetics  and  dynam¬ 
ics  of  large  molecules  excited  in  the  energy  range  0.2-6 
eV. 1-4  Photoselective  molecular  excitation  in  jets  above  6  eV 
is  expected  to  open  new  horizons  in  the  areas  of  extravalence 
spectroscopy,  nonreactive  and  reactive  intramolecular  dy¬ 
namics.  and  state  selective  molecular  photochemistry.5  A 
significant  step  in  this  direction  was  undertaken  recently  by 
three-  and  four-photon  ionization  studies6-9  in  jets  and  by 
femtosecond  time-resolved  investigations10  of  Rydberg 
states.  We  wish  to  report  on  the  merger  between  the  tech¬ 
niques  of  vacuum  ultraviolet  (VUV)  spectroscopy  and  super¬ 
sonic  jets,  which  allows  for  the  interrogation  of  absorption 
spectra  of  large  molecules  in  supersonic  expansions  in  the 
near  VUV  (6-9  eV)  region.  The  techniques  of  absorption 
spectroscopy  in  seeded,  planar,  pulsed,  supersonic  expan¬ 
sions"-1  ’  were  adopted  for  VUV  absorption  studies.  Planar 
supersonic  jets  of  Ar  seeded  with  benzene,  benzene-D6  or 
naphthalene  were  expanded  from  a  nozzle  slit  (dimensions 
0.27x90  mm,  repetition  rate  9  Hz  and  gas  pulse  duration 
300  /zsl  and  the  characteristic  rotational  temperatures  were 
Tr  ~20K.  Light  (spectral  range  1500-2000  A)  from  a  modi¬ 
fied  short-arc  high-pressure  Xe  lamp  (EG  &  G  Model  FX 
193U),  which  was  fitted  with  a  sapphire  or  CaF,  window, 
was  focused  with  an/  =  4  CaF,  lens  onto  the  jet  parallel  to 
the  nozzle  slit  at  a  distance  x  =  8  mm.  The  light  beam  was 
split  by  a  CaF:  plate  and  monitored  by  two  photomultipliers 
coated  with  a  (tetraphenyl  butadiene)  convector,  both  before 
and  after  crossing  the  planar  jet.  The  absorption  AI  /I0  was 
normalized  to  the  incident  light  intensity. 

VUV  absorption  spectra  of  jet-cooled  aromatic  mole¬ 
cules  provided  information  on  intramolecular  dynamics  of 
odd-panty  Rydberg  states  and  an  interference  between  Ryd¬ 
berg  and  intravalence  excitations.  The  electronic  origin  of 
the  3  Px,y['A  )  Rydberg  excitation  [3£.y(0|]  of  benzene  and 


of  benzene-Z>6  (Fig.  1)  is  superimposed  on  the  'A  lg  — *  '£  ,u 
intravalence  transition. u- 15  The  energy  of  3Piy(0)  (55  877 
cm- 1  for  C6H6  and  55  946  cm- 1  for  C6D6)  is  in  accord  with 
one-1415  and  three-photon16  bulb  spectroscopy.  The  line 
shape  of  3/^(0)  is  a  Lorentzian  (Fig.  1)  with  a  width 
(FWHM)  of  A  =  32.7  ±3.1  cm-'  for  CbHh  and 
A  =29.0  +  2.0  cm-1  for  C6D5.  The  homogeneous  line- 
widths  £  and  lifetimes  r  obtained  after  correction  for  instru¬ 
mental  resolution  are  £  =  28  +  3  cm  - 1  (r  =  0. 19  +  0.02  ps) 
for  C6H5,  and  £=24  +  2  cm-1  (r  =  0.22  +  0.02ps)  for 
C6D5.  Our  value pf  £  for  iPxy  (0)  is  close  to  the  linewidths  of 
27-33  cm  - 1  reported6  for  the  vibronically  induced  two-pho- 
ton  transitions  of  this  Rydberg  transition.  A  small  ( ~  10%) 
normal  deuterium  isotope  effect  on  the  relaxation  rate  may 
exist,  which  is,  however,  within  our  experimental  uncertain¬ 
ty,  implying  that  the  electronic  relaxation  process  of  3 Pxy  (0) 
is  characterized  by  moderate  energy  gap(s).  The  homogen¬ 
eous  line  broadening  of  the  ‘A  lg  —  'Elu  intravalence  excita¬ 
tion  of  jet-cooled  benzene  exceeds  £  for  3  P  (0)  by  about  one 
to  two  orders  of  magnitude,  which  is  in  accord  with  old  argu¬ 
ments, 17  so  that  the  intramolecular  relaxation  of  the  intrava¬ 
lence  'E  excitation  presumably  occurs  within  ~  10  fs,  i.e., 
on  a  time  scale  comparable  with  a  vibrational  period.  The 
Lorentzian  line  shape  of  'A  ]g  — *  3 Px  y  (0)  provides  compell¬ 
ing  evidence  for  the  absence  of  Rydberg-intravalence  inter¬ 
ference  effects'51819  in  benzene,  which  is  in  accord  with  Rice 
et  al. 15  The  absence  of  Fano  interference  effects70  can  be 
rationalized  within  the  framework  of  a  random  coupling 
model31’"  for  the  iPxy{0)-'E  lu  interstate  coupling  terms 
and  for  the  lA  lg  transition  moments  to  the  homo¬ 

geneously  broadened  background  manifold,  which  erode  all 
interference  effects. 3  Pronounced  Rydberg-intravalence 
interference  effects'8  are  exhibited  in  the  spectrum  of  jet- 
cooled  naphthalene  (Fig.  2).  The  excessive  line  broadening  of 
the  resonance,  e.g.,  £  =80  cm-l(r=70  fs)  for  the 
R  5  +  1360  transition  (Fig.  2),  is  close  to  the  corresponding 


FIG.  I.  Absorption  spectra  of  the 
'A  —  iP,,[A  j„ )  electronic  ongin  of 
benzene  and  benzene-0,,  in  pulsed  planar 
supersonic  expansions  of  Ar  at  p  =  105 
Torr  The  spectral  resolution  is  0.45  A. 
The  base  line  is  marked  as  zero  on  the  in¬ 
tensity  scale.  The  solid  curves  represent 
fits  of  the  line  shape  L  U  I  to  a  Lorentzian 
superimposed  on  a  constant  background 
L  |4  |  =  A  +  C/[U  -  a„|;  +  (4  /2|:).  For 
C„H,..  .4  =  5.7,  C=I7.2  (in  arbitrary 
unitsl,  A„=  1789.65  A  and  4  =  1.08  A 
(33.7  cm  ’ ')  while  for  C„D,,:  A  =  6.5, 
C=  1 2.6  (in  arbitrary  unitsl,  A, ,  =  1787.45 
A  and  4  =  0.88  A  127.5  cm  '). 
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FIG.  2  The  absorption  spectra  of  naphthalene  in  pulsed  planar  expansions 
of.Arip  =  105  Torn  in  the  region  1630-1650  A.  The  spectral  resolution  is 
0  5  A.  The  antiresonances  corresponding  to  the  R  ’  -r  1360  and  to  the  R  ‘ 
Rydberg  states  (according  to  the  assignment  of  Refs.  15  and  19]  are  marked 
by  arrows. 

room  temperature  value  r  =15  cm-1,15-19  whereupon  the 
excessive  homogeneous  broadening  make  the  thermal  con¬ 
gestion  effects.  Our  spectroscopic  approach  can  readily  be 
extended  throughout  the  LiF  region  (up  to  1 1  eV)  and  will 
contribute  towards  the  elucidation  of  the  microscopic 
aspects  of  intramolecular  dynamics  and  photochemistry. 
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Recent  measurements  of  energy-resolved  photoisomerization  rates  of  jet-cooled  diphenylbutadiene  molecules  are  analyzed  in 
terms  of  RRKM  calculations.  Perfect  agreement  over  a  wide  energy  range  is  obtained.  A  comparison  of  the  thermally  averaged 
rates  with  liquid-phase  data  is  made. 


1 .  Introduction 

The  photoisomerization  reactions  of  trans-stilbene 
and  all-trans-diphenylbutadiene  in  recent  years  have 
been  intensively  studied  under  isolated-molecule  con¬ 
ditions  in  low-pressure  gases  and  in  supersonic  beams, 
and  in  the  liquid  phase.  In  the  gas  phase,  one  of  the 
important  aspects  was  the  question  whether  statistical 
unimolecular  rate  theory  can  quantitatively  account 
for  the  energy-resolved  rates  of  the  non-radiative  pro¬ 
cesses  involved.  In  the  liquid  phase,  basic  theories  for 
barrier -overcoming  in  viscous  media,  such  as  Kramers’ 
theory,  were  tested.  In  order  to  conduct  a  meaningful 
analysis  of  the  liquid -phase  data,  a  direct  comparison 
with  the  gas-phase  data  appears  obligatory.  In  Kramers’ 
theory  [1 J ,  the  low -viscosity  limit  of  a  double-well 
isomerization  reaction  is  given  by  the  conventional 
“limiting  high-pressure  rate  constant”  of  a  thermal 
unimolecular  reaction 

/  mk(E)dE ,  (1) 

Eo 

with  the  thermal  equilibrium  distribution /(E),  the 
threshold  energy  Eq  ,  and  the  specific  rate  constant 
k{E).  Energy-resolved  gas-phase  experiments  give  di¬ 
rect  assess  to  k(E).  In  Kramers’  theory,  viscosity  ef¬ 


fects  in  the  liquid  produce  a  decrease  of  the  thermal- 
ized  rate  constants  below  the  maximum  value  of  eq. 
(1).  It  is,  therefore,  an  important  question  whether 
thermal  averaging  of  the  specific  rate  constants  k(E) 
from  isolated-molecule  experiments  provides  an  upper 
limit  for  the  liquid-phase  rate  constants. 

This  question  was  addressed  in  a  previous  article 
by  one  of  us  [2]  on  the  photoisomerization  of  trans- 
stilbene.  By  combining  the  isolated-molecule-energy- 
resolved  photoisomerization  results  from  jet-cooled 
[3—5]  and  thermal-distribution  [6,7]  experiments,  it 
was  possible  to  obtain  an  optimized  RRKM  fit  which 
perfectly  represented  the  experimental  data  offc(£). 
Calculating  k„  for  this  RRKM  fit  led  to  the  surprising 
result  that  liquid-phase  data  were  larger,  and  not  smal¬ 
ler,  than  by  up  to  one  order  of  magnitude.  Pos¬ 
sible  explanations  for  this  phenomenon  are  [2] ,  e.g., 
solvent  shifts  of  the  0—0  transition,  the  threshold  ener¬ 
gy  Eq,  or  the  activated  complex  frequencies,  which 
even  for  small  shifts  have  a  large  influence  in  low-bar¬ 
rier  unimolecular  reactions,  or,  alternatively,  incom¬ 
plete  vibrational  relaxation  of  the  lowest  vibrational 
levels  in  the  liquid  phase. 

Obviously,  a  similar  analysis  is  of  interest  for  the 
photoisomerization  of  diphenylbutadiene.  There  are 
recent  energy-resolved  measurements  of  k(E)  in  jet- 
cooled  molecules  by  Zewail  and  co-workers  [8]  over 
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the  range  0—1800  cm-1  and  by  Amirav  et  al.  [9]  over 
the  range  0-7500  cm-1 .  On  the  other  hand,  there  are 
low-viscosity  liquid-phase  results  by  Fleming  and  co- 
workers  [10,11]  over  a  wide  temperature  range.  A  pre¬ 
liminary  comparison  of  gas-  and  liquid-phase  data  in 
ref.  [12]  was  conducted  only  in  a  qualitative  way.  In 
our  present  work,  at  first,  we  try  to  obtain  an  optimized 
RRKM  fit  of*(£)  over  the  full  energy  range  0-7500 
cm-1  of  the  gas-phase  experiments.  With  these  results, 
we  calculate  k^( T)  and  compare  with  the  experiments 
in  liquid  ethane  by  Courtney  and  Fleming  [  1 1  ],  in  or¬ 
der  to  compare  with  the  gas-phase/liquid-phase  discrep¬ 
ancies  realized  in  the  trans-stilbene  system. 

2.  Specific  rate  constants  k  (E) 

We  do  our  RRKM  fit  for  the  gas-phase  experiments 
in  the  same  way  as  described  in  ref.  [2]  for  trans-stil¬ 
bene.  We  use  the  complete  set  of  84  excited-state  mo¬ 
lecular  frequencies  calculated  by  Pierce  and  Birge  [13] . 
For  simplicity,  we  fit  the  activated-complex  frequen¬ 
cies  by  multiplying  all  37  frequencies  related  with 
butadienic  carbon  atoms  (such  as  indicated  in  table  2 
of  ref.  [13])  with  a  common  scaling  factor  F.  We  use 
mode  84  (57  cm-1)  as  being  closest  to  the  reaction 
coordinate.  As  observed  before  [2,14] ,  scaling  of  a 
set  of  frequencies  or  of  a  single  frequency  does  not 
result  in  different  fits  ofk(£).  (In  trans-stilbene  either 
1 5  ethylenic  carbon  related  frequencies  or  a  single  fre¬ 
quency  out  of  the  72  molecular  oscillators  could  be 
scaled  [2] .)  The  threshold  energy  was  varied  around 
the  apparent  threshold  of  1050  cm-1  which  was 
“visually”  observed  in  ref.  [8],  Eq  and  the  scaling  fac¬ 
tor  F  were  the  only  two  parameters  varied  to  obtain  a 
fit  with  the  measured  k(E )  curve.  Since  there  were  no 
internal  rotors,  the  Beyer-Swinehart  counting  algo¬ 
rithm  with  an  energy  graining  of  1  cm-1  provided  an 
exact  count  of  activated-complex  number  of  chan¬ 
nels  W{E  -  Eq)  and  molecular  density  of  states  p(E) 
entering  into  the  standard  statistical  formula 

k(E)  =  W{E  -  E0)lhp(E ) .  (2) 

For  this  molecule,  W(E  -  Eq)  contains  a  factor  of  2 
for  the  symmetry  number. 

Table  1  and  fig.  1  give  k(E)  calculations  for  various 
parameter  sets.  An  optimum  agreement  is  obtained 
with  Eq  =  1 100  cm-1  and  F  =  1.11.  F  values  larger 


Table  1 

RRKM  calculations  for  the  photoisomerization  of  diphenyl- 
butadiene:  specific  rate  constants  k  (£)  and  high-pressure  rate 
constants  k^{T)  (for  molecular  and  activated-complex  pa¬ 
rameters,  see  text) 


£(  cm"1) 

*(£)(  s’1) 

£  (cm  *) 

k  (£) (s~‘ ) 

1100 

5.01  X  107 

1600 

1.42  x  109 

1150 

1.10  X  10® 

2000 

3.60  x  109 

1200 

1.64  x  10® 

2500 

7.06  x  109 

1250 

2.65  X  10® 

3000 

1.10  X  10‘° 

1300 

3.67  X  10® 

3500 

1.50  X  1010 

1350 

5.06  X  10® 

4000 

1.91  X  1010 

1400 

6.70  X  10® 

5000 

2.71  X  1010 

1450 

8.28  X  10® 

6000 

3.46  X  1010 

1500 

1.00  X  109 

7000 

4.14  X  1010 

1550 

1.21  X  109 

8000 

4.77  X  1010 

Jt«  (260  K) 

=  3.27  X  109  s"1 

k„  (290  K) 

=  5.80  X  109  s" 

k»  (270  K) 
k»  (280  K) 

=  4.02  X  109  s'1 
=  4.86  X  109  s"1 

k»  (300  K) 

=  6.83  X  109  s'1 

0  OX)  2000  1000  6000  8000 


E  /cm'1 

Fig.  1.  Specific  rate  constants  k  (£)  for  isomerization  of  ex¬ 
cited  diphenylbutadiene  (supersonic  beam  experiments  from 
ref.  [8]  (o)  and  ref.  [9]  (•);  optimized  RRKM  fit  from  this 
work  (full  line, £0  =  1 100  cm"1 ,  £  =  1.11),  with  the  calcu¬ 
lated  threshold  rate  constant  X). 
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than  unity  indicate  a  "more  rigid”  set  of  activated- 
complex  frequencies  as  compared  with  the  molecular 
frequencies.  A  similar  extent  of  rigidity  as  in  the  pres¬ 
ent  case  was  observed  in  trans-stilbene  where  F  -  1 .2 
was  fitted  with  about  half  the  number  of  ethylenic 
carbon  related  molecular  modes  [2],  Fig.  1  demon¬ 
strates  nearly  perfect  agreement  between  statistical 
unimolecular  rate  theory  (here  in  rigid  RRKM  form) 
and  experiment.  We  emphasize  again  that  this  remark¬ 
able  agreement  does  not  necessarily  imply  complete 
intramolecular  randomization,  since,  in  principle, 
minor  non-random  effects  could  be  compensated  by 
the  parameter-fitting  procedure.  Until  complete  infor¬ 
mation  concerning  the  potential  surfaces  of  the  sys¬ 
tem  becomes  available,  the  RRKM  analysis  cannot 
provide  a  conclusive  answer  to  the  cardinal  issue  of 
randomization.  Nevertheless,  the  optimized  RRKM 
calculation  provides  a  quantitative  fit  of  the  data. 
Since  the  fit  is  done  by  the  variation  of  two  param¬ 
eters,  E0  and  F,  an  independent  fit  of  E0  is  not  possi¬ 
ble.  However,  a  value  of 


E0  =  1 100  ±  50  cm-1 


For  this  reason,  it  appears  to  be  preferred  to  rely  on 
the  threshold  energy  Eq  derived  from  a  fit  of  the  full 
k(E )  curve  as  given  by  eq.  (3)  rather  than  on  the  “vi¬ 
sual”  threshold  energy. 


3 .  Comparison  of  gas-  and  liquid -phase  data 


probably  encompasses  the  widest  uncertainties  possi¬ 
ble. 

Whereas  the  data  are  perfectly  represented  over  the 
energy  range  1150-7500  cm-1 ,  there  is  a  minor  prob¬ 
lem  at  energies  below  1150  cm-1 .  The  detailed  mea¬ 
surements  of  quantum  yields  and  non-radiative  rates 
in  this  range  [8,9]  indicate  a  break  in  the  k(E)  curve 
near  1050  cm-1  with  k^^E)  s*  0  below  Eq  =  1100 
cm-1 .  At  present  it  is  difficult  to  provide  an  unambi¬ 
guous  interpretation  of  the  low-energy  tail  of  the 
knr(E)  curve.  The  marked  break  at  1050  cm-1  indi¬ 
cates  the  presence  of  an  additional,  slow  parallel  non- 
radiative  channel,  e.g,,  intersystem  crossing,  competing 
with  the  isomerization.  On  the  other  hand,  there  might 
be  rotational  effects  in  the  k(E,J)  curve  near  the  thresh¬ 
old  energy.  Such  effects  have  been  explicitly  demon¬ 
strated,  even  for  a  rigid  RRKM  case,  in  the  dissociation 
H2CO  -*  H2  +  CO  [15].  The  effects  are  most  pronounced 
for  small  /  near  the  threshold  energy.  Since  the  residual 
/-distribution  in  the  supersonic  beam  experiments  is  not 
exactly  known,  it  is  difficult  to  provide  a  quantitative 
account  similar  to  ref.  [15].  However, /-effects  and  the 
presence  of  a  second  channel,  probably  an  intersystem 
crossing  [10],  may  account  for  the  minor  discrepancies 
with  the  k(E)  calculation  near  the  threshold  energy. 


In  previous  liquid -phase  work,  Courtney  and 
Fleming  [11]  attempted  to  search  for  the  low -viscos¬ 
ity  turnover  of  the  Kramers’  theory  [1]  by  measuring 
the  photoisomerization  rate  in  liquid  ethane  and  pro¬ 
pane  at  temperatures  between  260  and  297  K.  In  these 
experiments  the  viscosity  and  the  temperature  were 
simultaneously  varied  so  that  viscosity  effects  and 
changes  in  the  Boltzmann  factor  are  hard  to  disen¬ 
tangle.  Measurements  with  varying  viscosity  for  one 
solvent  at  constant  temperature,  such  as  performed  by 
Brey  et  al.  for  stilbene  [16]  in  high-pressure  liquids,  or 
by  Troe  and  co-workers  [17,18]  with  halogens  in  sol¬ 
vents  with  densities  varying  continuously  from  the 
low-pressure  gas  phase  to  high-pressure  liquid  condi¬ 
tions,  are  not  available  for  diphenylbutadiene.  There¬ 
fore,  an  unambiguous  analysis  of  the  liquid-phase  data 
alone  appears  impossible. 

Our  present  calculation  of  specific  rate  constants 
k(E)  permits  us  to  construct  the  corresponding  high- 
pressure  rate  constants  k„  via  eq.  (1).  Table  1  includes 
k"(r)  values  for  our  optimized  RRKM  model.  A  com¬ 
parison  with  the  experimental  data  of  Courtney  and 
Fleming  [  1 1  ]  in  liquid  ethane  is  made  in  fig.  2. 


3.5  3.7 

tOO  K/T 


Fig.  2.  Comparison  of  “high-pressure  rate  constant”  k (this 
work)  with  measurements  of  thermalized  photoisomerization 
rate  constants  fcm  for  diphenylbutadiene  in  liquid  ethane 
(experimental  points  •  from  ref.  [11]). 
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Quite  unlike  trans-stilbene,  where  the  non-radiative 
rate  constant,  km,  in  liquid  solvents  was  found  to  be 
markedly  larger  than  k „  (up  to  one  order  of  magni¬ 
tude),  for  diphenylbutadiene  kni,  in  liquid  C2H6,  is 
slightly  below  .  It  also  shows  a  markedly  different 
temperature  coefficient.  Whereas  is  represented 
(between  260  and  300  K)  by 

=  8.3  X  10u  exp(-1000  cm-1  hcIkT)  s-1  .  (4) 

km  in  liquid  ethane  (between  263  and  291  K)  was 
found  [11]  to  be  equal  to 

km  =  5.8  X  1013  exp(-1854  cm'1  hc/kT)  s~l  .  (5) 

The  marked  differences  in  apparent  activation  energies 
suggest  that  km  for  liquid  ethane  is  still  in  a  range  of 
marked  viscosity  effects  far  below  the  low-viscosity 
maximum.  It  seems  to  us  that  the  low-viscosity  maxi¬ 
mum  of  fcnr  is  higher  than  k„  also  in  diphenylbuta¬ 
diene  indicating  solvent  shifts  in  the  molecular  param¬ 
eters  of  excited  diphenylbutadiene.  Only  experiments 
in  the  gas  phase  can  further  elucidate  this  point.  Ethane 
is  a  particularly  suitable  solvent  for  such  studies,  since 
density  can  conveniently  be  changed  in  a  continuous 
manner  from  the  low-pressure  gas  to  liquid  conditions 
[17].  Such  experiments  are  underway  [19] . 

4.  Conclusions 

We  have  shown  that  an  optimized  RRK.M  fit  pro¬ 
vides  an  excellent  representation  of  energy -resolved 
photoisomerization  rates  of  excited  diphenylbuta¬ 
diene  in  jet-cooled  free  molecules.  Thermal  averaging 
of  the  specific  rate  constants  k  (£)  leads  to  the  high- 
pressure  limiting  rate  constants  k„ .  A  comparison  of 
k„  with  liquid-phase  data  shows  markedly  different 
activation  energies.  There  are  strong  indications  that, 
as  in  trans-stilbene,  solvent  shifts  of  excited-state  mo¬ 
lecular  parameters  occur.  Before  an  analysis  of  liquid- 
phase  results  in  terms,  for  example,  of  Kramers’  the¬ 
ory  can  be  made,  these  solvent  shifts  have  to  be  de¬ 
termined.  Measurements  in  one  solvent  over  large  den¬ 
sity  ranges,  from  low-pressure  gases  up  to  high-pres¬ 
sure  liquids,  have  to  be  conducted  before  a  more 
meaningful  analysis  of  liquid-phase  data  appears  pos¬ 
sible. 
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The  quantum  path-integral  molecular-dynamics  method  was  applied  to  explore  the  structure,  en¬ 
ergetics,  and  dynamics  of  an  excess  electron  interacting  with  an  alkali-halide  cluster.  Four  distinct 
modes  of  electron  localization  were  established,  which  depend  on  the  cluster  composition,  size,  and 
structure;  they  involve  an  internal  E-center  defect,  an  external  surface  state,  dissociative  detach¬ 
ment  of  an  alkali  atom,  and  structural  isomerization  induced  by  electron  attachment. 


PACS  numbers:  71  45.Nt.  36.40  -t-d,  61  20.Ja 
Structural,  electronic,  dynamic,  and  chemical 
characteristics  of  materials  depend  primarily  on  the 
state  (phase)  and  the  degree  (size)  of  aggregation. 
Small  clusters,  i.e.,  finite  aggregates  containing  3-500 
particles,  exhibit  unique  physical  and  chemical 
phenomena,  which  are  of  both  fundamental  and  tech¬ 
nological  significance,  and  provide  ways  and  means  to 
explore  the  "transition”  from  molecular  to  con¬ 
densed-matter  systems.1  Theoretical  studies  of  clus¬ 
ters  were  hampered  by  the  relatively  large  number  of 
particles,  which  renders  the  adaptation  of  molecular 
science  techniques  rather  cumbersome,  while  the  lack 
of  translational  symmetry  inhibits  the  employment  of 
solid-state  methodology.  Molecular-dynamics  (MD) 
simulations,  consisting  of  the  generation  of  phase- 
space  trajectories  via  the  numerical  integration  of  the 
(classical)  equations  of  motion  for  a  many-particle  sys¬ 
tem,  are  particularly  suitable  for  the  study  of  the  struc¬ 
ture  and  dynamics  of  small  clusters  2,3  In  this  context, 
localized  excess  electron  states  in  clusters4  are  of  con¬ 
siderable  interest  with  regard  to  the  (nonreactive  and 
reactive)  mechanisms  of  electron  attachment,  the  for¬ 
mation  of  bulk  or  su'face  states,  and  the  role  of  the 
excess  electron  as  a  probe  for  the  interrogation  of  the 
nuclear  dynamics  of  the  cluster.  Futhermore,  quan¬ 
tum  phenomena  are  expected  to  be  pronounced  in 
such  systems  since  the  electron  wavelength  is  compar¬ 
able  to  the  cluster  size.  In  this  paper  we  report  on  the 
structure,  energetics,  and  dynamics  of  alkali-halide 
clusters  (AHC)  studied  with  classical  MD,  and  of  elec¬ 
tron  alkali  clusters  studied  with  the  quantum  path- 
integral  MD  method  (QUPID).5-8  AHC  were  chosen 


since  the  nature  of  the  interionic  interactions  is  well 
understood  and  in  view  of  the  abundance  of  model  cal¬ 
culations  and  experimental  information  of  these  sys¬ 
tems.9  Our  QUPID  calculations  establish  four  modes 
of  localization  of  an  excess  electron  in  AHC;  (i)  an 
E-center  defect  with  the  excess  electron  replacing  an 
internal  halide  ion;  (ii)  a  new  surface  state,  i.e.,  a 
“surface  E  center”  of  the  excess  electron;  (iii)  dissoci¬ 
ative  electron  attachment  to  AHC  resulting  in  the  for¬ 
mation  of  an  “isolated”  alkali  atom;  (iv)  structural 
isomerization  induced  by  electron  attachment.  Our 
calculations  establish  the  compositional,  structural, 
and  size  dependence  of  these  various  localization 
mechanisms. 

The  QUPID  method6-8  was  applied  to  a  system  of  an 
electron  interacting  with  an  AHC  consisting  of  A  ions 
(N,  cations  and  iV2  anions).  The  interionic  potential 
energy  within  the  AHC  is  FAHC- '£'u(t>ij(Rij),  with 
the  interionic  pair  potentials  4>u(/?w)  being  given  by 
the  Born-Mayer  potential  with  the  parameters  deter¬ 
mined  by  Fumi  and  Tosi.10  The  electron- AHC  poten¬ 
tial  is  ^(rj-S/d’f/tr-R/)  consisting  of  a  sum  of 
electron-ion  potentials,  which  are  described  by  the 
purely  repulsive  pseudopotential  <t>fl(r)  -  e2/ r  for  the 
electron-anion  interaction  and  by  the  local  pseudopo¬ 
tential11  <P,/(  r)  =  -e2/Rc,  r  ^  Rc,  and  4>,/(r) 
-  -  e2/ r,  r  >  Rc<  for  the  electron-positive-ion  in¬ 
teraction.  The  Hamiltonian  is  //-  Ke  +  V,  -I-  KAHC 
+  Fahc>  where  K,  and  KaH c  are  the  kinetic-energy 
operators  for  the  electron  (mass  m)  and  of  the  ions 
(masses  Mt=M ,  and  A/2),  respectively.  Observables 
are  obtained  from  the  quantum  partition  function 
Z  -  lim,_  .[Z^with 


Z,  -  Tr(exp(  -  r  K  t,HC)exp(  -  r  K,  )exp(  -  t  KAHC)exp(  -  t  V, )  ], 


where  r  -  /3/p and  /3  —  ( kBT)  - 1  is  the  inverse  temperature.  If  we  make  use  of  the  free  particle  propagator12  Zp  is 


z,-  n 

*■1,2 
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where 


Mr  =  X  ^^y'r,-r(  +  |)-+ L,(r,)/p 


N  "  dM , 

Mr  =  X  X  T-3-j-(R/ui-R/(/  +  i))2+  Vahc/p- 

/_[,■_!  m  p 

Equation  (1)  maps  the  quantum  problem  onto  the 
classical  statistical  mechanics  of  .V  +  1  particles,  each 
consisting  of  a  periodic  chain  (necklace)  of  p  pseu¬ 
doparticles  (beads)  with  nearest-neighbor  harmonic  in¬ 
terchain  interactions,  whose  strengths  depend  upon 
the  masses  ( m ,  Mu  A/2>,  the  temperature  (T),  and 
the  pseudoparticle  number  (p).  When  the  thermal 
wavelength  [kT- (fift2/ Mr)]/2]  is  smaller  than  any 
relevant  length  scale,  the  Gaussian  factor  in  Zp 
reduces  to  a  delta  function  and  the  necklace  collapses 


to  a  classical  particle.  This  is  the  case  for  the  ionic  part 
of  our  system.  The  average  energy  of  the  system  is 
given  at  equilibrium  by 

P 

E  =  4^-  +  <  ^AHC)  +^f  +  P_1(X  Mr/)) 

W  /-I 

with  the  electron  kinetic  energy 

K  =  3  1  r/W  (r  _r  A 

K'  2/3  2  p~x\  8r,  '  p)! 

which  consists  of  the  free-particle  term  (3/2/3)  and  the 
contribution  from  the  interaction  (A'.m)  with  the 
ions.13  The  statistical  averages  indicated  by  angular 
brackets  are  over  the  Boltzmann  distributions  as  de¬ 
fined  in  Eq.  (1).  This  formalism  is  converted  into  a 

numerical  algorithm  by  noting14  the  equivalence  of  the 
sampling  described  above  to  that  over  phase-space  tra¬ 
jectories  generated  via  MD  by  the  classical  Hamiltoni¬ 
an 


[  pm(Ti-Tl  +  i)2 

Mr ,) 

l  2fr2*2 

P 

i- i  1  i- l  ^  /-i  m 

the  mass  m‘  being  arbitrary  and  taken  as  m'=  1  u. 

Numerical  simulations  were  performed  for  an  elec¬ 
tron  interacting  with  sodium-chloride  clusters  at  about 
room  temperature.  On  the  basis  of  examination  of  the 
stability  of  the  variance  of  the  kinetic-energy  contribu¬ 
tion  £im,  the  number  of  “electron  beads”  was  taken 
as  p  =  399.  By  use  of  an  integration  step  of 
Ar  =  1.03x  10“ 15  sec,  long  equilibration  runs  were  per¬ 
formed  [(l-2)x  1 04 A  f ] .  The  reported  results  were 
obtained  via  averaging  over  8x  103Ar,  following  equili¬ 
bration.  The  electron-ion  pseudopotential  parameters 
were  varied  by  changing  the  cutoff  radius  Rc  in  the 
range  (3.22-5.29)a0-  From  QUPID  calculations  on  a 
single  Na  atom,  the  atomic  ionization  potential  is 
reproduced  (see  Table  I)  for  Rc  =  5.29 a0,  which  seems 
to  be  a  too-high  value  for  the  characterization  of  this 
pseudopotential.6  Therefore,  on  the  basis  of  pseudo¬ 
potential  parametrization  studies11  and  a  recent 
QUPID  study8  of  F centers  in  molten  and  solid  KC1  in 
which  the  same  form  of  pseudopotential  was  em¬ 
ployed,  a  value  of  Rc  =  3.22a0  (yielding  a  value  of 
-0.3005  a.u.  for  the  electron  binding  energy  to  Na+, 
see  Table  I)  is  preferred.  Fortunately,  our  conclusions 
regarding  internal  and  surface-localization  modes 
remain  unchanged  with  respect  to  reasonable  varia¬ 
tions  of  this  parameter.  In  simulations  involving  the 
electron-AHC  interaction,  different  initial  conditions 
were  employed,  two  of  which  for  the  lNa14Cl13  J  ■*"  sys¬ 
tem  are  portrayed  in  Figs.  1(a)  and  Kb).  The  result¬ 
ing  final  state  of  the  system  was  found  to  be  indepen¬ 
dent  of  the  initial  conditions. 

It  has  previously  been  suggested  on  the  basis  of 
zero-temperature  structural  calculations,9  and  has  been 
confirmed  by  our  classical  MD  simulations,  that  when 


the  size  of  the  clusters  increases  (;V3s20),  the  NaCl 
cyrstallographic  arrangement  is  preferred  for  particular 
stability  for  clusters  forming  rectangular  structures 
even  if  the  number  of  positive  and  negative  ions  is  not 
equal.  Therefore,  we  have  chosen  to  study  first  the  in¬ 
teraction  of  an  electron  with  [NanC^]'1'  and 
[Nai4Cl12]  +  +  clusters,  which  exhibit  pronounced  sta¬ 
bility.  In  Figs.  1(c)  and  1(e),  we  present  our  results 
(using  Rc  =  3.22a0)  for  the  equilibrium  electron- 
charge  distribution  obtained  from  2D  projections  of 
the  necklace  edge  points,  and  for  the  nuclear  con¬ 
figuration  of  the  clusters.  In  both  cases  the  electron, 
which  starts  in  either  initial  configuration  as  shown  in 
Figs.  1(a),  Kb),  has  been  localized.  However,  two 
distinct  modes  of  electron  localizations  are  exhibited 
involving  internal  and  external  localization  for  the  dou¬ 
bly  charged  and  singly  charged  cluster,  respective¬ 
ly.  The  vacancy-containing  configuration  of  the 
[Nai4C!!2]  +  +  cluster  stabilizes  an  internally  localized 
excess  electron  state,  with  the  e  surrounded  by  six 
Na+  ions  in  an  octahedral  configuration  and  by  twelve 
Cl"  ions  [Fig.  1  (e) ]  which  is  similar  to  the  case  of  an  F 
center  in  the  extended  solid.  The  electron  affinity  of 
the  cluster  EA  =  £J  +  Ec  is  obtained  by  summing  ihe 
electron  binding  energy 

£'!=^  +  *,n,  +  />-1i<F<(r,)> 

1- 1 

and  the  cluster  reorganization  energy  Ec=(VAHC) 
-  (  FAHC)0,  where  (  Fahc>o  is  the  potential  energy  of 
the  "bare”  AHC  in  the  absence  of  the  electron.  The 
ionic  configuration  of  the  e-[Na!4CI!2]+  +  cluster  is 
somewhat  distorted;  however,  the  gain  in  £|(- 0.249 


Volume  54.  Number  16 


PHYSICAL  REVIEW  LETTERS 


22  Apr  11  14.S5 


TABLE  I.  Average  equilibrium  temperature  ((D),  interionic  cluster  potential  energy  ((UAHC)).  electron  interaction 
kinetic  energy  (£',„,).  electron  kinetic  energy  A',  =  3/2/3+  A,nt,  e-AHC  interaction  potential  energy  ((  (',)),  electron  binding 
energy  (£J),  cluster  reorganization  energy  (£f),  electron  affinity  of  cluster  ( EA  ),  and  Cartesian  components  of  the  “electron 
necklace”  gyration  radii  ( Rj.Ry.R,7)  Atomic  units  are  used  (energy  and  (  T)  in  Hartrees,  length  in  Bohr  radii).  Variances  are 
given  in  parentheses.  Calculated  values  of  {T)  and  (  Lahc)  for  the  “bare”  clusters:  [NauCluL  0.976(0.1 1 1) x  io~3 
-3.7296;  [NauCI,j]  +  .  0.948(0. 106)  x  lO’3,  -3.5911;  [Nal4Cl12]+ +,  0.954(0.093 )  x  10'3,  -3.3151;  [Na5CI4]  +  .’ 

1.01  (0.38) x  10-3,  -1.0999. 


e  -l  Na14CI|3  J 
Rc  —  3.22 
e-[Na14Cl,j] 
Rc  =  5.29 
e-[Na14Cli2] 
Rc  =  3.22 
e-[NauCln] 
Rc  =  5.29 
e-[Na5CI4]  + 
Rc  =  3.22 
e-[Na5Cl4]  + 
Rc  =  4.36 
e-Na+ 

Rc  =  3.22 


Rc  —  4.36 
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Re  -  5.29 
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.983  -3.4856 

.035) 
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.031) 
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(0.956) 

2.6389  0.0278  -0.3012  -0.2734  0.0940 

(0.582) 
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FIG.  1.  Ionic  configurations  and  “electron  necklace"  distributions  for  an  excess  electron  interacting  with  sodium-chloride 
clusters.  Small  and  large  spheres  correspond  to  Na+  and  Cl-  ions,  respectively.  Dots  represent  2D  projections  of  the  “elec¬ 
tron  beads.”  (a),(b)  Alternate  initial  configurations  of  lNat4CI|3  j  *  +  e  [in  (b)  the  e  bead  is  localized  to  the  right  of  the  cluster] 
which  achieve  the  equilibrium  configurations  corresponding  to  surface  states  given  in  <c),  Rc-3.22o0,  and  (d)  /?c-  5  29a0 
(e)  Equilibrium  configuration  of  e-[Na,4Cl12]  *  +  .  Rc  =  3  22o0.  exhibiting  an  internal  £  center,  (f)  Initial  state  of  <’-[NasCl4] +  . 
(g)  Equilibrium  configuration  of  [Na?Cl4 ] +  with  Rc  -  3  22a0  resulting  in  dissociation  of  Na,  (h)  Equilibrium  configuration  of 
e-(NajCI4]*  with  Rc  =  4.36a0  which  corresponds  to  structural  isomerization. 


au.)  exceeds  the  loss  in  Ec  ( 0.0799  a.u.)  favoring 
internal  localization  (Table  I).  It  is  of  interest  to  note 
that  the  total  energy  of  ( Na!4Cl |  > ]  +  +  is  rather  close 
to  that  of  lNaMCI|il  +  ,  so  that  the  electron  binding  en¬ 
ergy  in  the  cluster  is  similar  to  that  of  a  negative  ion, 
analogously  to  the  situation  for  F-center  formation  in 
(extended)  ionic  crystals,8  thus  establishing  the  domi¬ 
nance  of  short-range  attractive  interactions  in  electron 
trapping  (localization)  phenomena.  A  drastically  dif¬ 
ferent  localization  mode  is  obtained  in  the  <>- 
[Na14Clu]+  system  (Fig.  1  (c)  1 ,  where  a  surface  state 
is  exhibited.  For  Rc  =  3.22a0  the  electron  localizes 
around  an  Na+  surface  ion  (Fig.  1(c)],  leaving  an 
essentially  neutral  (NanCI|j]  cluster,  which  interacts 
with  the  (partially)  neutralized  Na  atom  mainly  via  po¬ 
larization  of  the  electron  cloud  (with  a  residual  ionic 
binding  of  -0.002  a.u.).  We  refer  to  this  state  as  a 
cluster-surface  localized  state,  which  bears  close  analo¬ 
gy  with  Tamm’s  crystal-surface  states. IS  When  the 
pseudopotential  cutoff  Rc  increases,  the  surface  state 
becomes  more  extended  (see  Fig.  1  (d)  for  Rc 
=  5.29a0)  (the  choice  of  an  abnormally  high  value  for 
Rc  (i.e.,  5.29a0)  prevents  internal  localization  in  the 
[Na|4Cl12l  +  +  cluster  since  the  internal  region  is  then 
predominantly  repulsive  (see  Table  I),  resulting  in  a 
surface  state]  A  measure  of  the  spatial  extent  of  the 
localized  electron  is  given  by  the  gyration  radius  of  the 
“electron  necklace"  R2  -  ( \/2p2)  <  v(,(r, •- tj)2) . 
which  demonstrates  (Table  I)  the  enhanced  localiza¬ 
tion  in  the  e-(Na|4CI|2]+  +  system  and  the  anisotropy 
of  the  electron  distribution  in  the  e- ( Nat4CI ! 3 1 +  sys¬ 
tem.  An  estimation  of  the  extent  of  the  electron  ther¬ 
mal  wave  packet  is  obtained  also  by 

p  1/2 

-f-j -  2<<r*-ri-i)2> 

P  1  i- 1 

which  for  a  free  electron  at  room  temperature  yields 
Rt-J 3Xr—  56 a0  (A^—  32.34o0).  In  all  the  calcula¬ 
tions  reported  herein  RT  values  of  52a0-55a0  were  ob¬ 
tained,  i.e.,  the  same  as  the  free-electron  value  to 
within  statistical  significance  (compare  to  typical  in- 
lerionic  distance  of  5a0  in  NaCl  clusters). 

In  smaller  clusters  novel  effects  of  dissociative  elec¬ 
tron  attachment  and  cluster  isomerization  induced  by 
electron  localization  will  be  manifested.  We  have 
studied  the  smallest  singly  ionized  cluster  exhibiting 
high  stability,  i.e.,  (Na5CI4]  +  ,  for  which  the  lowest  en¬ 
ergy  configuration  is  planar  with  four  Na+  ions  at  the 
corners  and  one  at  the  center  of  an  approximate 
square.’  This  cluster  possesses  an  isomeric  structure 
(less  stable  by  0  014  a.u.)  in  which  the  ions  are  ar¬ 
ranged  in  a  distorted  pyramidal  configuration.’  Adding 


an  electron  to  the  ground-state  planar  configuration  of 
(Na5CI4)+  (with  the  Na+  pseudopotential  being 
characterized  by  Rc  =  3.22a0)  transforms  the  system 
into  a  neutral  (Na4CI4]  cluster  with  a  planar  ring  struc¬ 
ture  and  a  dissociated  neutral  Na  atom  (Table  I  and 
Figs.  1(f)  and  Kg),  corresponding  to  the  initial-  and 
final-state  configurations,  respectively].  To  demon¬ 
strate  that  this  process  is  driven  by  the  localization  of  e 
around  a  single  Na+  ion,  we  have  performed  further 
simulations  by  decreasing  the  cation-electron  binding 
orength,  taking  Re-4.26a0.  In  this  case,  extreme  lo¬ 
calization  is  not  sufficiently  counterbalanced  by  e  bind¬ 
ing.  Instead,  the  planar  structure  transforms  to  the 
isomeric  pyramidal  configuration  with  the  electron  lo¬ 
calized  as  a  diffuse  cloud  about  the  tip  of  the  pyramid 
(Fig.  1(h)].  Electron  localization  accompanied  by 
structural  isomerization  will  constitute  a  prevalent 
phenomenon  for  AHC  with  smaller  e-cation  binding 
energy,  i.e.,  the  heavier  alkali  metals.  In  view  of  the 
intimate  interrelationship  between  structural  isomeri¬ 
zation  and  melting  of  clusters,16  it  will  be  interesting  to 
explore  the  melting  of  such  finite  systems  induced  by 
electron  localization. 
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